arXiv: 1506.04253v2 [gr-qc] 16 Sep 2015 


Parameterized post-Newtonian limit of Horndeski’s gravity theory 


Manuel Hohman 
Teoreetilise Fuiisika Labor, Fiiusika Instituut, 

Tartu Ulikool, Ravila 14c, 50411 Tartu, Estonia 

We discuss the parameterized post-Newtonian (PPN) limit of Horndeski’s theory of gravity, 
also known under the name generalized G-inflation or G^-inflation, which is the most general 
scalar-tensor theory of gravity with at most second order field equations in four dimensions. 
We derive conditions on the action for the validity of the post-Newtonian limit. For the most 
general class of theories consistent with these conditions we calculate the PPN parameters 
7 (r) and P(r), which in general depend on the interaction distance r between the gravitating 
mass and the test mass. For a more restricted class of theories, in which the scalar field is 
massless, we calculate the full set of PPN parameters. It turns out that in this restricted 
case all parameters are constants and that the only parameters potentially deviating from 
observations are 7 and /3. We finally apply our results to a number of example theories, 
including galileons and different models of Higgs inflation. 



I. MOTIVATION 


The most striking observations in modern cosmology are the accelerating expansion of the 


universe 


background 



whose cause has been named dark energy, and the homogeneity of the microwave 
which is conventionally attributed to an inflationary expansion of the very early 


universe [Y|, |^J. However, both dark energy and inflation are yet unexplained phenomena. A large 
and important class of theories aiming to explain these phenomena is based on the introduction of 
a scalar field mediating gravity in addition to the usual metric degrees of freedom . Various 


theories belonging to this class have been successfully applied to cosmology [1^, [H 

An important feature of general relativity, which one wishes to retain also in scalar-tensor 
theories of gravity, is the fact that its gravitational field equations contain at most second order 
derivatives of the dynamical fields. This restriction is imposed since higher derivative theories will 
in general lead to instabilities and ghosts jl^ . One may therefore ask which is the most general 
class of scalar-tensor theories whose field equations are of at most second order. This question 
was answered already by Horndeski [l^, although his work did not receive much attention until 
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recently and was re-derived in a different, but equivalent formulation 


18l-l2ll|. Since then, many 


particular theories 
and dark energy 


aelonging to this class have been studied, in particular as models of inflation 


22 


However interesting its performance in cosmology, a viable gravitational theory must also pass 
the tests on local scales, e.g., give a good account of the motions in our solar system. A natural 
framework for such a check is the parameterized post-Newtonian (PPN) formalism 


23!, 


2|. It 


characterizes gravity theories by a set of ten parameters, which have been measured with high 


precision in various solar system experiments 


25l433l| . Through the availability of this high precision 


data, the PPN formalism has become an important testbed for the viability of gravity theories. 
This work extends and generalizes an earlier result on the PPN parameters 7 and /3 for a 


class of scalar-tensor theories of gravity with a general potential in the Jordan frame 
have also been calculated in the Einstein frame taking into account screening effects 
terms of invariants under conformal transformations and scalar held redehnitions 


341, which 
351 ] and in 


36l | . Also the 


solar system physics of other theories belonging to the Horndeski class of gravity theories have 
already been studied and it has been argued that the “hfth force” mediated by the scalar degree 
of freedom should be suppressed in order to reproduce the observed general relativity limit. For 


this purpose, several screening mechanisms have been studied 2^, such as the chameleon 37h39|, 
symmetron [40l. l41 1 or Vainshtein 42h44| mechanisms, and it has been shown that these mechanisms 


can achieve consistency of the theory with solar system observations 


4 ^|. 


In this work we complement these studies by an analysis of theories in which screening mech¬ 
anisms do not play a significant role, so that the standard PPN formalism can be applied. We 
explicitly calculate the PPN parameters 7 and j3 for a general class of Horndeski theories, and 
the full set of PPN parameters for a more restricted class with a massless scalar field, in order to 
show that also in this case consistency with solar system observations can be achieved, without 
employing any screening mechanisms. 

The outline of this article is as follows. In section m we display the action and discuss the 
structure of the field equations. In section IHII we expand these field equations in a weak field 
limit around a Minkowski background. The post-Newtonian limit of this expansion is discussed in 
section lYl The post-Newtonian gravitational field equations are then solved for a static point mass 
source in section El which yields the PPN parameters 7 and /3. The full set of PPN parameters is 
obtained in the case of a massless scalar field in section IVll In section I VIII we compare this general 
result to current observations of the PPN parameters. We apply our findings to a few example 
theories in section IVHIl and end with a conclusion in section IIXI 
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II. ACTION AND FIELD EQUATIONS 


In this section we provide a brief overview of the structure of the action and the field equations 
of Horndeski’s gravity theory. The starting point of our derivation is the action, which takes the 
form [l^ 

5 

'S' = 'y ^ / d g£,i[gfj_u^ 0] + *Sm; Xm] • (1) 

i=2 

Here Sm denotes the matter action and Xm collectively all matter fields. The gravitational part of 
the action, which depends on the metric and a single scalar field (f>, is given as an integral over 
the four-dimensional spacetime manifold, where the Lagrangian is composed of the terms 

C2 = K{4>,X), C3 = -G3{(j),X)D(j), C4 = Gi{(P,X)R + Gix{(l>,X)[{D(f)f, 

A = - ^G 5 x(</>, X) [{ncj)f - 3(n</))(V^V,<?i)2 + 2{V^V,ct>f] . (2) 

Here we introduced the notation 


(V^ V,</))3 = V^V,</.V" V VVa , A = - i (3) 


for the d’Alembert operator □ and derivatives of the scalar field, and indices are raised and lowered 
with the metric g^y. The functions K,G 3 ,G 4 ,G 5 are free functions of the scalar field (p and its 
kinetic term X. Each choice of these functions determines a distinct gravity theory. We denote 
derivatives of these functions by a subscript, e.g., G^x = dG^jdX. 

The gravitational field equations are derived from the action ([1]) by variation with respect to 
the metric and the scalar field. It follows from the structure of the action that the field equations 
take the general form 




( 4 ) 


i=2 


i=2 


i=2 


where is the energy-momentum tensor of the matter fields Xm- The terms 0^^, and are 
obtained from the variation of the different Lagrangians in the gravitational part of the action. 
Their full form is rather lengthy and listed in the appendix of [l^. However, for practical purposes 
it turns out to be easier to replace the metric field equation with its trace-reversed analogue 

5 


E’^: 

i=2 


flU 


-It -1 

— — r, 


^g^yT 


( 5 ) 
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where the trace-reversed metric terms are given by 

= (6) 

These are the field equations we will be working with in this article. For the purpose of calculating 
their post-Newtonian limit, we first need to bring them into a more manageable form using a 
perturbative expansion around a fixed background solution. This will be done in the next section. 


III. PERTURBATIVE EXPANSION 


In order to calculate the parameterized post-Newtonian limit of Horndeski’s gravity theory we 
will need a perturbative expansion of the field equations, which we displayed in the preceding 
section, around a fixed background solution. This background solution will be given by a fiat 
Minkowski metric and a constant cosmological background value of the scalar field, so that 
the perturbative expansion assumes the form 

9tiu = + hij,,y, 0 = 4>-h'(/;, X =. (7) 


Besides assuming that the background is homogeneous and isotropic, we thus also assume that 
it is stationary, i.e., constant in time. The physical reasoning behind this assumption is that 
we particularly consider the situation at or close to a fixed point of the background evolution of 
the scalar field, so that we can neglect any effects from a dynamical background. For the post- 
Newtonian limit it will be necessary to expand the terms (and thus also and up 

to the quadratic order in the perturbations and V' around this background. This will be done 
in this section. 

Recall that the action of Horndeski’s gravity theory, and thus also the field equations, depends 
on the choice of four free functions K, G 3 , G 4 , G 5 , which depend on the scalar field (p and its kinetic 
term X. Their Taylor expansion around the cosmological background value $ takes the form 

CXD 

Ki<P,X)= ( 8 ) 

77T,,n=0 


where the coefficients R(m,n) are given by 

^ Qm+n 


K, 


(m,n) 


m\n\ dcp'^dX'' 


-K{4>,X) 


(9) 


0=4.,x=o 


and similarly for the remaining functions G 3 ,G 4 ,G 5 . Each term is of the order 

Q^^m+ 2 n-^^ Using these expansions, the terms constituting the field equations listed in Q take 
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the form 


^/ii/ ~ 2'^(2i0)^a‘^^ 

+ ^^( 0 ,i)^/^;^ 5 pV’ 5 '’V’ - > ( 10 a) 

Gftu = G?,{i,o)d^i'4^dy'ilj - ]^G^^ifl)r]^ydpi:dPi: , ( 10 b) 

^fiu = G4(0,0)G^i/[li^] + G'4(l,0)^p!^nV’ — G4 (i^o)^M^i/V’ + G'4(i^o)G'Ati^[l^^]V’ + G4(0,0 )G^!^[/i^] 

+ 2G4(2,0)^P!^V'I^V’ + G^4(1,0) h^i> ni, + G4(i,o)rv[/i']5pV' - G4(i,o)hM.^""r"p,.[/^']5rV’ 

- 2G4(2,o)V’5p'9,,V' - G^i^Q^rjpi.hp^dPd'^'ip - G4(o,i)nV’'9;,5i.V’ + G4(o,i)«9p5^V’^^5i.V^ 


+ 2^4(2,o)??pi/5pV'5^V' + ^G4 (o,i)??p^ [(□V’)^ - dpdcr'ipdPd'^'il;] - 2 G^ 2 ,o)df^'>pduip , (10c) 

= G5(i,o)nV'5p5i/V' - G^i^Q)dpdpil}dPdyil} - ^G5(^o)^pi^ , (lOd) 

for 

= -^( 1 , 0 ) + ‘^^{2,0)4’ + 3-?^(3,o)V’^ “ {1,1)9p'G'd^'G , (11a) 

= 2^3(2,o)9pV’'9^V’, (11b) 

= G4(i^o)-^[^^] + 2^4(2,o)-R[^^]V’ + G4 (i^o)-^[ 1^^] + b^4(l,l) [(OV’)^ “ SpSo-V’^^^^^V’] ) (Hc) 

p5 = o, (lid) 

for P^ and 

Jp = -K{o, 1 ) 9^14 - ^(i,i)V’^p'0 , (12a) 

Jp = 2G3(i_o)5pV' + 4^3(2,o)V’^pV' + ^^3(0,1) [2nV’5p'0 - dpidpiidPip)] , (12b) 

= 2 G 4 (o,i)Gp.[h']a> - G4 (i,i) [ 2 Dil 2 dpiP - dpidpipd^iP)] , (12c) 

j5 = -2G5(i,o)Gp,[hi]a>, (12d) 


for Jp up to the quadratic order in hpi, and ill. Note that here we have changed our notation 
from the one we used in the previous section. From this section to the end of section lYD where 
we discuss perturbations around a flat background, □ = rj^^'^dpdi, denotes the flat Minkowski 
d’Alembert operator and indices are raised and lowered using the flat metric rjp^. Further, we have 
introduced the notation F[h^] for the term which is of order n in the expansion of F with respect 
to the metric perturbation hp^. For our calculation we further need the trace-reversed terms W'py 
dehned in equation (l 6 |) and the divergences V^Jp. From the expansions (|10p and (|12p one easily 
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derives the expansions 




2^i0,0)'ntJ.iy + 2 -^( 0 , 0 )^!^ + + -K(^iQ^h^j,v'4) + 


2K{o,i)dfM'ipdu'ip 
T^lv = Gzi^dpii’dy'il), 

74 


(13 


la 


(13b) 


2^4(1,— ^4(1^0) *9/^91/V’ + G^4(l,0)-^/i!^[^^]V’ + G'4(o,0)-R/i!^[^^] 

“ <^4(2,0)— -G4(i^O)^/i!^l^'0 + 2^4 


— <^4(0,0)“ 2^' 


2^4(1,0)ViMuhpad^d'^i; 

- ‘^G^2,0)^dpd^'il; + G4(i^0)r^Mi^[^^]^pV' - G^Q^i^a^pdpd^'ip + G^o^i^dpdp^pdPd^'ip 

- G4(^2,o)Vt^udp'iljdP'ilj - 20^(2^)9^11^0^11) + ]^Gi(i^Q)rip^r]P^T'^p„[h^]dril) , 




nL = G. 




for and 

J2 = -K( 0 ,i)n^ + K( 0 , 1 ) + K(0,i)r?^"r%,5pV’ - i^d.DV’QV^ - i^(i,i)9pV^5^V’, (14a) 

v^^Jl = 2G3(i,0)nv^ - 2G3 (i, 0) - 2G3(i,0)??^"r%,apV’ + 4G3(2,0)V^nV’ 

+ 4^3(2,+ G'aio,!) [{Giipf - dpdaipd/'d^il)] , (14b) 

= 2G^(o,i)GpAh^]d^d''il) - 2G^^^^^ - dpd^il)df>d^il)] , (14c) 

V^4 = -2G5(i,o)Gp,[/r']a^9'^V’. (14d) 


for V^ J)j. From these expressions we see that the only Taylor coefficients relevant for our discussion 
will be 


-^(0,0) ’ -^(1,0) ) -^(2,0) ’ -^(3,0) ’ -^(0,1) > -^(1,1) ’ ^3(1,0) ’ G^3(2,0) > 

^3(0,1) ) ^4(0,0) ’ ^4(1,0) > G'4(2,0) > ^4(0,1) > ^4(1,1) ) ^^5(1,0) ' (1^) 

All other terms in the Taylor expansion would lead to corrections of at least cubic order in the 
field perturbations. 

The perturbative expansions of the field equations around a fixed background spacetime dis¬ 
played in this section will be a central ingredient for the calculation presented in this article. The 
second main ingredient will be the parameterized post-Newtonian formalism, which is built upon 
these perturbative expansions and will be discussed in the following section. 
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IV. POST-NEWTONIAN APPROXIMATION 

The main tool we use in this article is the parameterized post-Newtonian (PPN) formalism 
, which we briefly review in this section in the context of the given scalar-tensor theory of gravity. 
The key idea of the PPN formalism is the assumption that the matter which acts as the source 
of the gravitational field is given by a perfect fluid, whose velocity in a particular, fixed frame of 
reference is small, measured in units of the speed of light, and that all physical quantities relevant 
for the solution of the gravitational field equations can be expanded in orders of this velocity. We 
will now show how this expansion in velocity orders proceeds for the quantities we need in our 
calculation in the following sections. 

The starting point of our calculation is the energy-momentum tensor of a perfect fluid with rest 
energy density p, specific internal energy 11 , pressure p and four-velocity u^, which takes the form 

= (p + pn + p)uf^u'^ + . (16) 

The four-velocity is normalized by the metric so that u^u'^g^i, = — 1. We will now expand 
all dynamical quantities in orders 0{n) oc |i;|"' of the velocity u* = v}/vP of the source matter in a 
given frame of reference, starting with the field variables. For the metric g^^ this is an expansion 
around a flat Minkowski background, 

= 9^11' + + 0{5 ), (17) 

while the scalar held 4> is expanded around its cosmological background value, 

(/) = $ + '!/) = ^> + + 0(5). (18) 

Here each term resp. is of order 0(n). In order to describe the motion of test bodies 

in the lowest post-Newtonian approximation an expansion up to the fourth velocity order 0(4) 
is sufficient. A detailed analysis shows that not all components of the metric and the scalar held 
need to be expanded to the fourth velocity order, while others vanish due to Newtonian energy 
conservation or time reversal symmetry. The only relevant, non-vanishing components of the held 
variables are given by 

^00 > > ^00 > ■ (19) 

In order to determine these components for a given matter source we must assign velocity orders 
also to the rest mass density, specific internal energy and pressure of the perfect fluid. Based on 
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their orders of magnitude in the solar system one assigns velocity orders 0(2) to p and IT and 0(4) 
to p. The energy-momentum tensor ()16p can then be expanded in the form 

21)0 = p ^1 + n -|- + 0(6), (20a) 

Toj = -pvj + 0(5), (20b) 

Tij = pviVj + p6ij + 0(6). (20c) 


For later use we also expand the trace-reversed energy momentum tensor introduced in the field 
equations ([S|) into terms of velocity orders and obtain the expressions 

Too = \p+ \p^ + - \phoo ’ (21a) 

Toj = -pvj + 0(5), (21b) 

Tij = ^pSij + pviVj -h ^Ph^if - ^P^ij + 0(6). (21c) 


We further assume that the gravitational field is quasi-static, so that changes are only induced by 
the motion of the source matter. Time derivatives do of the metric components and other fields 
are therefore weighted with an additional velocity order 0(1). 

In order to solve the gravitational field equations, which inherit a gauge symmetry from the 
diffeomorphism invariance of the gravitational action, we finally need to fix a gauge for the metric 
tensor. A useful choice for the class of scalar-tensor theories we consider can be constructed in 


analogy to the gauge condition introduced in 


48| and takes the form 


1 


1 . 


G. 


4(1,0) 


o ^jj,i + r, ^00,i — ^ 

^ ^ '-"4(0,0) 


1 


G, 


4(1,0) 


' 0 , 1 ) hojj hjjfl — 

^ <-"4(0,0) 


■ 0 , 


0 • 


In this gauge the Ricci tensor up to the required order takes the form 


( 22 ) 


f) __ 

-n-OO — 2 00,kk o 00,kk^ 

Ga 


G. 


4(b0) ,( 2 ) 


n '",00 
<-"4(0,0) 


<^ 4 ( 1 , 0 ) ( 2 ) , ( 2 ) _ 1 ^( 2 ) ^( 2 ) 

+ 2 ^ 4(0 0 ) 2 oo,j%o,j + 2 ’ 

1 . 1 . io\ G, 


p _ _ ^/,(^) I ’~^4(i,o) iX2) pj(K\ 

-flOj — 2^0j,kk ^%0,0j T ^,0i ’ 


Ri. = 


Ga 


2 d,kk ' 


4(0,0) 


(23a) 


(23b) 

(23c) 


We now have all necessary tools at hand in order to derive the post-Newtonian limit of Horn- 
deski’s gravity theory. Using the expansions (El) and (|18l) we can derive the post-Newtonian 
gravitational field equations from the expansions shown in section m keeping only the terms (|19p . 
We further apply the chosen gauge by inserting the gauge fixed Ricci tensor (|23l) into the resulting 
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equations. Finally, we can solve these equations following increasing velocity orders. This will be 


done first under the assumption of a static, spherically symmetric gravitationa 
a single point mass in the following section, following the approach detailed in 


field generated by 


34l | . For a more re¬ 


stricted class of theories, corresponding to a massless scalar field, we can go beyond this assumption 
and fully solve the post-Newtonian gravitational field eqnations, and will do so in section IVll 

V. STATIC SPHERICALLY SYMMETRIC SOLUTION 

We will now construct a solution to the post-Newtonian gravitational field equations derived 
from the perturbative expansions detailed in the previous sections. The central assumption in this 
section will be that the source of the gravitational field is given by a single point mass, so that the 
gravitational field is static and spherically symmetric. The general form of this field is explained in 
section IV A 1 The solntion will yield three physical quantities which we will subsequently determine: 
the effective gravitational constant Geg in section IV Bl and the PPN parameters 7 in section IV Cl 
and j3 in section IV D1 The calculation presented here is carried out in full analogy to the calculation 
displayed in an earlier work for a more restricted class of scalar-tensor theories fl- 

A. General form of the solution 


The starting point of our calculation is the assumption that the source of the gravitational field 
is given by a single point-like mass M, whose energy-momentum tensor is of the form (I16p with 


p = M6{x), n = 0, p = 0, Vi = 0. 


(24) 


This simple matter source induces a static and spherically symmetric metric, which can most easily 
be expressed using isotropic spherical coordinates. In the rest frame of the gravitating mass we 
use the ansatz 


goo = -1 + 2Geff(r)U(r) - 2Glf,{r)P{r)U\r) + + 0(6), 

90j=O{5), 

Qij = [1 -f 2Ges{r)-f{r)U{r)] dij + 0(4). 


(25a) 

(25b) 

(25c) 


where r denotes the radial coordinate and the spherically symmetric, static Newtonian potential 
U{r) is given by 

M 


U{r) = 


(26) 
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In the potential we collect terms of order 0(4) which are not of the form such 

as the gravitational self-energy. Note that we have changed the normalization of the Newtonian 


potential (I26p in comparison to previous work j^ . where we used an additional factor Here 


we have omitted this factor in order to be consistent with the standard normalization given in 
which will be relevant in section lYD 


231, 


The metric (I25p contains three unknown functions which we need to determine. These are the 
effective gravitational constant Geff(r) and the PPN parameters 7 (r) and /3(r). The latter two can 
be defined either as the coefficients of the effective gravitational potential t/efr = as shown 

in the metric (l2^ or as the coefficients 7 eff = GefT 7 and /3efr = of the Newtonian potential 

terms U and t/^. The first definition invokes the interpretation that the measured values of 7 and 
j3 can be related to the effective gravitational potential Ues, while the second definition suggests 
to relate the measured values of 7 efr and /3efr to the Newtonian potential t/ of a fixed mass M. 
We choose the hrst definition in this article since the mass of the Sun, which dominates the solar 
system physics, is determined from its gravitational effects on the planetary motions. 

The Newtonian potential U (r) we introduced here is of second velocity order, so that the zeroth 
velocity order solution is given by the flat Minkowski background = 7 ^ 1 /. For this solution to 
be consistent with the gravitational field equations, we must check that it solves the zeroth order 
in their perturbative expansion, which corresponds to the equations of motion for the cosmological 
background. The corresponding equations follow from our assumption that this background is 
stationary and are given by 


2^(o,o)VtJ.iy — 0 ) -^( 1 , 0 ) — 0 • 


(27) 


As one can see, these are solved only if the Taylor series coefficients 77 ( 0 , 0 ) -^( 1 , 0 ) vanish. We 

will therefore restrict ourselves to theories which satisfy these conditions. 


B. Newtonian approximation 

We will now determine the effective gravitational constant Gefr(r) appearing in the metric 
solution . The starting point of this calculation is given by the gravitational field equations (jH) 
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and ([5]). At the second velocity order they are given by 

= G4(o,o)-Rto + 2^4(1,) (28a) 

\lf = 04(0,0)df - ^04(1,„,v,<2% - G4(i,„,V,f , (28b) 

0 = (-^( 0 , 1 ) “ 2G3(i^o)) + 2 i^( 2 ,o)V’^^^ + ^4(1^0)-^*'^^ ) (28c) 

( 2 ) 

where Riiv is the second velocity order part of the Ricci tensor. Here we already used the restriction 
R(o,o) = -^( 1 , 0 ) = 0 derived from the zeroth order field equations. The components of the trace- 
reversed energy-momentum tensor are given by 

T^o=\p^ = (29) 

In order to eliminate the second order Ricci scalar = —R^q+R^^'^ from the third equation (I28cp . 
we can take the trace 


= 04(0, 0 ) (-4? + df) - 304(,,„,«,2 . 

over the first two equations (I28all and (I28bp . This yields the scalar field equation 


-^( 0 , 1 ) “ 2G3(i^o) + 3 


^4(1,0) 




G, 


( 2 , 0 )' 


4(1,0) 


^ ^*4(0,0) ^ 

which has the form of a screened Poisson equation, 

'^fi - = -Ci,p 

where we introduced the constants 

G. 


2G. 


4(0,0) 


p = 0, 


= 


\ 


-2K, 


(2,0) 


K, 


( 0 , 1 ) 


- 2G 


3(1,0) 


1^2 

+ 37^ 


C-Tp — 




2G. 


K, 


’G 


4(0,0) 


( 0 , 1 ) 


-2G 


3(1,0) 


+ 3 


fi2 

4(1,0) 


G. 


4(0,0) 


(30) 


(31) 


(32) 


(33) 


4(0,0) 


In order for a consistent solution to exist which is compatible with the perturbation ansatz we 
restrict ourselves to the case A"( 2 ,o) ^ 0 O'!!!! -^(o,i) “ 2 G 3 (i_o) + o)/^4(o,o) > 9- The solution 

is then given by 


47rr 


(34) 


for the point mass source 

41^—^—1. f 2 ^ 

In the next step we use equation (|28al) and insert the second order Ricci tensor Rqq from 

I ( 2 ^ ( 2 ^ j—i 

equation (I23al) . In the resulting equation for we eliminate the term using equation (l32|) 
and finally obtain 


hm,ii = - C2P 


(35) 
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where the constants ci ^2 are given by 


Cl = -2 


C2 = 


G'4(l,0)-f^(2,0) 


G. 


4(0,0) 


-f^(o,i) “ 2G3 (i 0 ) + 3 


^4(1,0) 


a 


4(0,0) 


G. 


4(0,0) 


1 ^ ^4(1,0) 


2G, 


K, 


4(0,0) 


( 0 , 1 ) 


-2G 


3(1,0) 


+ 3 


^4(1,0) 


-1 


G. 


4(0,0) 


(36a) 

(36b) 


For the point mass source (1241) this equation has the solution 




C2 + ^(e-”^^'--l) 


mz 


(37) 


By comparison with the metric component (|25a|) we read off the effective gravitational constant 


Geff(r) - 


C2 + 


*"l*"b („—m^r _ 


m 


-ie- 


1 ) 


'll) 


(38) 


In order to interpret this result for G^s as an effective gravitational constant we need to choose an 
experiment in which the gravitational interaction takes place at a constant scale r = tq. We can 
then choose units in which Geff(ro) = 1. This corresponds to a rescaling of the parameter functions 
K, Gs, G 4 , G 5 . However, we cannot make this choice globally, and hence cannot remove the factor 
Gefr(r) from the metric (1251) b v a choice of units in which Gefi = 1 , as it is conventionally done in 


the basic PPN formalism 


23|. This is the reason for the ambiguity in the definition of the PPN 


parameters 7 and (3 we discussed above. 


C. PPN parameter 7 (r) 


We now come to the calculation of the PPN parameter 7 , which can be read off from the 

(‘2'\ 

spatial perturbation component h\j , as it appears in the corresponding term in the metric (I25cl) . 
For this purpose we use the field equation (I28bl) and insert the second order Ricci tensor Rlj 
from equation (|23cp . As we did in the previous section when calculating the effective gravitational 
constant, we eliminate the term using equation (I32p and finally obtain 

- C4p) Sij , (39) 


where the constants C 37 are given by 


C3 = 2 


C4 = 


G4(l,0)-f^(2,0) 


G, 


4(0,0) 


-^(0,1) “ 2G3(i^o) + 3 


r-f2 

^4(1,0) 


G, 


4(0,0) 


G, 


4(0,0) 


^4(1,0) 


2G, 


K, 


4(0,0) 


( 0 , 1 ) 


-2G 


3(1,0) 


+ 3 


^4(1,0) 


-1 


G. 


4(0,0) 


(40a) 

(40b) 
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This equation has the solution 



M 

47rr 


C4 + ^(e-^ 


mz 


- 1 ) 


(41) 


for the point mass source (1241) . By comparison with the metric component (I25cp one reads off the 
PPN parameter 


7(r) = 


C2 + ^ (e-rnv-'- - 1) 2a; + 3 + 


(42) 


where the constant a; is given by 


UJ = 


^*4(0,0) 

9^2 

4(1,0) 


(-^( 0 , 1 ) ~ 2 ^ 3(1 0 )) • 


(43) 


The result reproduces a previously derived result for the PPN parameter 7 (r) for a more restricted 


class of scalar-tensor theories of gravity 


34l |. It thus also yields analogous limiting cases, which are 


obtained as follows. In the limit —)• 0 and fixed finite u, the PPN parameter 7 approaches the 

known value 


O' -|- 1 
a; -|- 2 


(44) 


for scalar-tensor gravity with a massless scalar field 
value 7 = 1, independent of m^. The same value 7 
a massive scalar field with m^r 1 . 


49|. In the limit uj ^ 00 we find the limiting 


= 1 is also approached in the limiting case of 


D. PPN parameter /3(r) 


We finally come to the calculation of the PPN parameter /3, which is read off from the component 
/igg^, which follows from the metric term (I25al) . Since the field equations at the third velocity order 


1^(3) 

2 


^4(0,0) 


R 


(3) 

Oi 




(45) 
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are solved identically for the static, spherically symmetric solution we consider here, we can directly 
proceed with the fourth order field equations. These take the form 


2^00* — C4(o,|))i?oo^ + 2^4(1,OjV'if - - 2^4(1,0)’/',00 + ^4(1,0)^00^'/'*^* 

+ G4(2,0)«g’V’<"> - io4,i,o)4oVf - io4,4,„,/!SfV';g’ - 
+ G4„ - iG4„ o, ( m - h,?> + hS’d .*?> . 


uu ,22 2 U >U 2 

+ - ^Q(r,o) + \h^l) ^ , 

- G'4(i,0)V',''^j^ + 2^(2,0) 

^<^3(1^0) - 2'^(°d) “ 2^4(2,0)^ ^ 


(46a) 


^ij 


0 


^ . *4(2,0) j + 2‘^4(l,0)V’,TO'^d' + ^4(1, 

- ^ij - 2^4(l,0)^lf V’,fci + 2^iil,0)h^kli’^kUij - 2^4(2, 

+ 5G4(i,o, (Mi)+Mh - Ci) MM + 5^4(1,0) (4h - MM*. 

- G'4(2,0)V',fcV,fc^<^ii + (<^5(1,0) - <^4(0,1)) (4’[fj'>P% - i’^fk 
= (K(0,1) - 2G3(1,0)) V’J^ + 2iM.o + G.n n^4R(") + 3K, 




i-V. 


,( 2 ) 


,jk) ’ 


(46b) 


+ 3iM(3^o) ^ 


+ 27^(2,0)V’^‘^^ + ^4(1,0)^^'^^ . - y^,^, y 

V’,S + 

+ -2G3(2.o)) (2V^gV^'^ +V’,?V,® 

(3G4(i,i)-G3(0,1)) 


+ (2G3(i^0) “ ^(0,1)) 

+ ‘ 


,( 2 ) 


+ 2 (<-^5(1,0) “ ^4(0,1)) ^if'^^ij ’ (^6c 


where the components of the trace-reversed energy-momentum tensor read 


'^00 = Ip^ + P^‘^ - \p^m + \p > + PViVj + \phf^ - ]^p6ij . 


(47) 


We can eliminate the fourth order Ricci scalar = —Rqq^ + R^m — Rif from 

the fourth order scalar equation (j46cp using a suitable linear combination of the fourth order 
equations (|46ap and (|46bh and the second order equations (|28ap and (|28bh . which reads 


^pU - = G4(o,0) (4'^ - 4? - 4^4? - 4o 4?) - 3^4(1,0)V^?? + 3G4(i,0)V^; 

+^^4{i,0)h[f^fj +‘2K(^2fi) 


.,( 2 ) 

,00 


+ 3^4(4 0) ( l ~ -4 I + -hn 


+ (^^3(1^0) “ 9.^(0,!) “ 6^4(2,0)^ ~ 6^4(2, 


+ ^4(1,0) (4? - 4?) + (^5(1,0) - G4(o,1)) - 44 ,^ 


.,( 2 ) 

4i 


( 48 ) 
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The resulting equation for the scalar field ■0^^^ finally takes the form 


V’fJ - + d2 ( 0 ^^^) + dslpf-l/jf + ^40^^^^^ + d5 ( 0 ® 

+ ^ 602 ' 0 ^^ + dr 0 ® + 

+ dgh^fijf^ + dioG^f+ diipU + dup , 


(49) 


V 


where the constants di,..., di 2 are listed in equation (lAip in appendix Terms of the forms 
/ ( 2 )\^ ( 2 ) ( 2 ) ( 2 ) ( 2 ) 

( 0 ii 1 ) 0V0V contain squared second derivatives of the Newtonian potential, and 

thus squares of the matter density. These terms do not appear in the standard PPN formalism 
and their influence on the current methods to measure 0 and other PPN parameters must be 
determined by a separate phenomenological discussion. We will not enter this discussion here, 
and therefore restrict ourselves to gravity theories in which the free functions in the action o are 
chosen so that d^ = d^ = dig = 0. By comparison with their values listed in equation ()A1[) this 
corresponds to the restrictions 


G^3(0,1) — 3G4(i^i) , (^4(0,1) — ^'5(1,0) 


(50) 


on the Taylor expansion coefficients of the functions G 3 , G 4 , G 5 . We now insert the point mass 
source (1241) and the already determined solution for 0^^0 and We further neglect terms of 
the form pU, which correspond to gravitational self-energies and thus contribute only to the term 
$(4)(r) 

in the metric component (|25al) . The resulting equation then reads 

^-m^r ^-2m^r ^-m^r ^-2m^r ^-m^r ^-2m^r 

ru - + 62^:^ + 63 ^;^ + 64 ^:^ + ’ ( 51 ) 


where the constants ei,... ,65 are listed in equation (IA2I1 in appendix O Prom this equation we 
obtain the solution 


0 (^)=/l^ + / 2 " 


-2m^r ^-m^r ^-m^r 

-h 03-In(m^r) -|- 04-Ei(-m.0r) 


r 

+ 05-Ei(- 2 m. 0 r) -h 06-Ei(-3m^r), 


(52) 


with constants fi,..., Jq listed in equation 
integral, which is defined by 


in appendix IS Here Ei denotes the exponential 


r 

Ei(-a;) = - / 

J X 



( 53 ) 
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In order to determine we now eliminate the term from equation (l46aD by making use of 
equation (|49D . The resulting equation then takes the form 

= 91^,S + 92 + 95 ) + 96^3^® 


/,(4) 


+ 97 (h\fi - V'!? + + 99/iS?V’S + 910^1?V'S 


(54) 


+ 9ll4o,i^!f + 912^0 V’,?? + 913 ^!)oV^^^ + 9i44o!i4Si + 915 /iSf 4o!ii 

+ + 9i7/5n + gispu^ + qi'iph^QQ + 920P, 

where the constants qi,..., (720 are listed in equation (IA4I) in appendix |A1 As it was also the case 

\ ^ (o^ (‘2^ (2^ (2^ 

in equation (@91), we find terms of the forms ['4’ii) j and Glj 'ip\j , which do not appear 

in the standard PPN formalism and which we therefore eliminate by the restriction 95 = 96 = 
qio = 0. A calculation of these coefficients shows that they already vanish as a consequence of the 
restriction ()50p we imposed earlier, so that all terms involving squares of second derivatives drop 
out. Into the remaining equation we insert the point mass (|24p and the previously determined 
solutions for the scalar field and the metric perturbations. Again we neglect all gravitational 
self-energy terms of the form pU . This yields us the equation 

3— 2m^r 


, ( 4 ) e e 

Kin = + S2 


—2m./,r 


1 

+ ^ 7—7 + ^8- 


r‘ 

^-m^r 


"ijjR g g 

+ sa —5 - 1 " S4“ 


—2m.i,r 


g g 

+ 'S 5 - 7 -1" S 6 — 


Q-m^r gin^r 

In(m^r) -|- sg-Ei(—m^r) + sio-Ei(—2m.0r) 


(55) 


+ Sii——Ei(-3m^r), 

where we used the constants si,..., sn listed in equation 
given by 


in appendix 13 The solution is then 


—2m^r 


, ( 4 ) 1 e e 

Ko =ui^ + U2 -n-h Us-^ 

1, 


T U 4 - 


^-m^r 


r 

Q-m^r 


+ U 5 - 


^-2m^r 


+ Ua- 


^-m^r 


In(m^r) 


-I- U7-Ei(—m^r) + U8Ei(—m.0r) + ug 


1, 


Ei(—m^r) + ttio-Ei(—2m.0r) (^6) 


+ niiEi(-2m^r) -|- ni 2 ^ Ei(-2m^r) -|- rtia——Ei(-3m^r), 

with constants rti,...,ui 3 listed in equation (|A6p in appendix [3 By comparison with equa¬ 
tion (|25ap and after inserting all coefficients listed in appendix 3 we finally read off 


/3(r) = I + 


I 


uj + T — 4uja _ 


(2UJ + 3 + e-™V''-)2 ] 2a; -b 3 


2m^r 


+ (2a; -b 3)m^r 


e In(m^r) — (m^r + Ei(—2m^r) — -e 


—2m-^jJr 


+ 


6pr -b 3(3a; -b r -b 6(T -b 3)m^r 
2(2a; -b 3)m^ 


[e"*’^^Ei(—3m^r) — e ”^’^'’Ei(—m^r)] 


( 57 ) 
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where oj is given by equation (|43p and we further introduced the abbreviations 
G^4(0,0)G^4(2,0) 


a = 


r'2 

4(1,0) 


r = 


'-'4(0,0) 


^*^ 4 ( 1 , 0 ) 


(-^(1,1) “ 4G'3(2,0)) > — 


^4(0,0)-^(3,0) 


^3 

^4(1,0) 


(58) 


We thus see that the result for j3{r) has essentially the same structure as a previously found result 


for a more restricted class of scalar-tensor theories of gravity 
behavior of the exponential integral in the case x S> 1 , 


It follows from the asymptotic 


/ 1! 2! 3! \ 

Ei(-x) 1-- + -- - + ... , (59) 

X \ X X-^ J 

that all terms involving cr, r or fall off proportional to and are thus subleading to the 

terms involving only uj and which fall off proportional to e”™'/''’. We therefore conclude that at 
large distances m^r )§> 1 from the source the contributions of a, r and ^ may be neglected. This 
means in particular that the comparison of j{r) and /3(r) with experiments in the large distance 
limit detailed in is valid also in the more general case of Horndeski’s gravity theory considered 
here. We will explain this limit in more detail in section IVIIl 

Again we consider the three limiting cases which we already discussed for 7 . In the limit 0 

and fixed finite uj we obtain 




^ ^ 00 + T — 4a;cr 

(2a;-k3)(2w + 4)2 ’ 


(60) 


which essentially reproduces the known result for a massless scalar field m . The second case 
Lo ^ 00 and arbitrary yields the limit /3 = I. We also find the limiting value /3 = I in the case 
fn^r I of a massive scalar field. 

This result completes our solution to the post-Newtonian field equations for a static point mass 
source of gravity. We have calculated the metric up to the first post-Newtonian order as displayed 
in equation (j25n . From our calculation we obtained expressions for the effective gravitational 
constant (j38p and the PPN parameters 7 (1421) and (3 (I57p . In the next section we will consider 
a more restricted class of theories, for which we can solve the post-Newtonian field equations for 
arbitrary mass distributions and obtain a full set of PPN parameters. 


VI. FULL SET OF PPN PARAMETERS FOR A MASSLESS SCALAR FIELD 

In the previous section we have determined the PPN parameters 7 and (3 from the static, 
spherically symmetric metric ansatz (1251) . For the consistency of this ansatz we had to impose the 
conditions (|^7P and ([^Up on the functions K,Gs,G 4 ,G 5 in the gravitational action. We have seen 
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that the parameters 7 and j3 depend on the distance between the mass source and the probing 
test mass, due to the fact that the scalar field acquires a non-vanishing mass (13,81) . We will now 
consider a further restricted class of theories in which this mass term vanishes, which is the case if 
the Taylor expansion coefficient K^ 2 fi) vanishes. Further, we require that the mass-like (derivative 
free) term in equations (m and (IMD vanishes, which is achieved by ^^(3,0) = 0 - When 

this restriction is imposed, it will turn out that the gravitational held equations can be solved 
for arbitrary matter sources given by the energy-momentum tensor (|16l) . and that their solution 


assumes the standard PPN form, from which the full set of ten PPN parameters can be read off 


23|. 


In order to determine this solution, we will solve the gravitational held equations by increasing 
velocity orders - the second velocity order in section IVI A[ the third velocity order in section IVIBI 
and the fourth velocity order in section IVI Cl 

A. Second velocity order 

We start by solving the gravitational held equations at the second velocity order following the 
same steps as in the preceding section, i.e., we hrst determine then /iqq^ and hnally For 
the scalar held we see that equation (f52|l reduces to 




(61) 


where is given by equation ()33p . and thus takes the form of an ordinary Poisson equation. The 
solution is given by 


= 2±u, 

dvr 


(62) 


where we have introduced the Newtonian potential 


U{t,x) = / d^x 


/ p{t, x') 


\x — X' 


(63) 


Note that for the point mass 


U reduces to the previously introduced spherically symmetric 

( 2 ) 


Newtonian potential (I2UI1 . Analogously, equation (I55|l governing reduces to 


J,( 2 ) 

Kin = -C2P , 


(64) 


where C 2 is given by equation (j36bl) . and thus has the solution 


, ( 2 ) _ 

Ko - 4 ^^- 


( 65 ) 
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We can compare this result to the corresponding metric component (I25ap in the spherically sym¬ 
metric case. From this we see that the effective Newtonian constant is given by 


G = 


Stt 


( 66 ) 


Here we have dropped the subscript “eff” in order to indicate that G is now really a constant, 
in contrast to being an effective quantity, which depends on the distance between the gravitating 
mass source and the test mass. 

We proceed by solving for the metric component . The corresponding equation (l39|l reduces 
to 


hfj,kk = -C4p6ij , (67) 

where C 4 is defined in equation (I40bll . This is again a Poisson equation, which is solved by 

4 ? = ^US,J . ( 68 ) 

With this result we have determined the metric at the second velocity order. 


B. Third velocity order 


(3) 

We now come to the third velocity order metric component %■, which is determined by equa¬ 
tion (145 p . In the case of a massless scalar field we consider in this section this equation reduces 
to 


'’Oljj Q 


1 C2 

-P'^i - . 


4(0,0) 


Stt 


(69) 


Note that in contrast to section FV PI this equation is not satisfied identically under the assumptions 

(3) 

made in this section, since we have not assumed that Hq-' and Vi vanish and U is time independent. 
We thus find that the solution is given by 


^ 

0 * - levr 


C2 


G. 


4(0,0) 


Hi - C2Wi 


(70) 


where we used the third order PPN potentials 

J \x-x'\ J 


,3 ,p{t, x')Vi{t, f') ^8 _ f ,3 - x'j) 


\X — X' 


•/|3 


(71) 


This result determines the metric at the third velocity order. 
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C. Fourth velocity order 


We finally come to the solution of the gravitational field equations at the fourth velocity order. 
This calculation is considerably simpler than the corresponding calculation in section IV Dl since 
from our restriction ^^( 2 , 0 ) = 0 follows that qiQ = 0, so that does not appear in the fourth order 
metric equation ([5T|) . We thus do not need to calculate and can directly proceed with solving 
for /iqq^ Inserting the solutions found for the second and third velocity order in sections IVI Al 
and IVI Bl equation (IMD reduces to 

^oo!m = wiUfio + W2U^iU^i + w^pU + Wipn + w^pv'^ + tcep, (72) 


where the constants wi,...,w^ take the values ()A7p listed in appendixThe solution is given by 




tC3 


104 


wq. 


Wl Wl 


2 dvr / \ dvr 

where the newly introduced PPN potentials are given by 


dvr 


dvr 


$ 


X — X 




J \x-x'\ J 


X 


X — X' 

,/ p{t, x') 


A{t,x) = J 


/m2 


/ [Vi{t ,x'){xi - x()] 


\x — X‘ 


•/|3 


, B{t,x)= [ d^x' {xi - x'i) 

J I TT* 'r' I 


X — X \ 

X — x'\ dt 


(74) 


With this result we have finally calculated all metric components to their respective velocity orders, 
which are required to determine the PPN parameters. 


D. PPN gauge and PPN parameters 


We now use the solution for the metric components /iqq\ h\^\ and calculated above 
in order to determine the PPN parameters. For this purpose, the metric must be in a particular 
gauge, in which it takes the form 231 


hf^ = 2U, (75a) 

= 2jU6,j , (75b) 

hgi = ~2 + 47 + q:i — a 2 + Cl ~ 2^)Vi — —(1 + 0(2 ~ Cl + 2^)Wi , (75c) 

= -2/31/2 - 2 Cd>u 7 + (2 + 27 + as + Cl - 2C)^i + 2(1 + 87 - 2/3 + C 2 + C )^2 ( 75 d) 

+ 2(1 + C3)$3 + 2(37 + 3C4 - 2C)$4 - (Cl - 20-4, 
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where in addition to the previously listed PPN potentials we have also included the Whitehead 
term 


$ 


w 




(rx"p{t, 


Xi 


\x — X‘ 


V|3 


Xj-Xj 

\x — x''\ 




(76) 


Note in particular that this metric must not contain the PPN potential B, in contrast to our 
result dZaD- This indicates that the solution we found is not yet in the PPN gauge. We thus 
need to eliminate the potential .S by a suitable gauge transformation, i.e., by a suitable change of 
coordinates. It turns out that the PPN gauge is achieved by introducing new coordinates x^ given 
by 


x^ = x^ + 


Wl 

?0 =- -^Xfi ; 


= 0 , 


where y is the superpotential 


xit,x) = — / d^x'p{t,^)\x — x'\ 


This gauge transformation changes the metric to 

Wl 

9oo — 9oo + + B — $i), Pol = goi + ~ 

In this new gauge we thus find the metric components 


) 9ij — 9ij 


^( 2 ) ^ S2ljj 

00 47r ’ 

m = ^us, 

u dvr 

= — 

Or 

W2 


ij , 


C2 - 


W5 


G. 


4(0,0) 


+ 


Wl 


Vi -{C2 + 


Wl 


hm = -(W2 + 


dvr 


dvr/ 


Wi 


\ ^ W4 We 

I d>9 ——d>3 - —<I>4 , 


dTT 


dvr 


(77) 

(78) 

(79) 

(80a) 

(80b) 

(80c) 

(80d) 


up to the required velocity orders. It is conventional to work in the normalization G = 1, in which 
one can directly read off the PPN parameters from the metric equation (1751) . This normalization 
is obtained by multiplying the gravitational part of the action and thus the free functions 
K,G 3 ,G 4 ,Ge, with the constant C2/87r. After applying this normalization we find the metric 
components 


Ag> = 2U , 

W = 2—U5., 






C2 

T 

2 


1] 5 


Wl \ 
G^4(0,0)C2 8C2/ 


+ 


= 327 r 2^[/2 - 2—$1 - 
Co C2 


Gdvr^ 


1 

2 8c2; 


]W,, 


W2 + 


Att J 


) $2 - 2—dia - 2—<I>4 . 


^W4 

C2 


C2 


(81a) 

(81b) 

(81c) 

(81d) 
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We can now compare this to the standard form (1751) of the PPN metric. Reading off the PPN 
parameters and inserting the previously introduced constants listed in appendixl^we finally obtain 


oj + 1 ^ w + T — Aoja 

a; + 2 ’ 4(w+ 2)2(20;+ 3) ’ 


«! = a2 = as = Cl = C 2 = Cs = C4 = ? = 0, 


(82) 


where we used the abbreviations (j43p and (j58p . This result shows that Horndeski’s theory of gravity 
belongs to the class of fully conservative theories, in which momentum and angular momentum are 
conserved, and in which there are no preferred-frame effects. Theories of this type are characterized 
by the PPN parameters ai = a2 = as = Ci = C2 = Cs = C4 = 0- Further, it is free of preferred- 
location effects, or Whitehead effects, which is indicated by the vanishing Whitehead parameter 
This leaves only the two PPN parameters 7 and /3 which potentially deviate from observations, as 
we will argue in the following section. 


VII. COMPARISON WITH OBSERVATIONS 

We now briefly compare the results obtained in the previous two sections to the values of the 
PPN parameters measured in solar system experiments. We restrict our discussion to the PPN 
parameters 7 and /3, since for a massive scalar field these are the only parameters we have calculated 
in section El while for a massless scalar field these are the only non-trivial PPN parameters, 
according to our calculation in section m where we have seen that all other parameters take the 
value 0 , in agreement with observations independently of the choice of a particular theory from 
Horndeski’s class. See Q for a recent review of the values of the full set of PPN parameters. 

We start our discussion with the case of a massive scalar field considered in section El Here 
we can essentially distinguish two regimes: a light scalar field / short interaction distance with 
m^r 1, and a heavy scalar field / long interaction distance with m^r S> 1. We will not consider 
the intermediate case, as it can simply be obtained by interpolation. In the limit of a light scalar 
field, the values of the PPN parameters approach their values in the massless case, which we will 
discuss later in this section, so that for now we will focus on the heavy scalar field case. Note that 
the PPN parameters 7 and /3 shown in (|42l) and (15711 depend exponentially on m^r. Keeping only 
the leading order terms we find 

^ = 1 - , /3 = l + ^p^ln{m^r)e-^^^+ ( 83 ) 

Zid + o ZlU + o 

so that the PPN parameters in this limiting case depend only on the constants and w. In 
order to derive bounds on these constants, we must consider measurements of 7 and /3 at a fixed 
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interaction distance r = tq. Currently the most stringent bounds of this type are obtained from the 
time delay of radar signals sent between Earth and the Cassini spacecraft on its way to Saturn 2^ . 


The experiment yielded the value 7 — 1 = (2.1 ± 2.3) -10 The radio signals were passing by the 
Sun at a distance of 1.6 solar radii or rg ~ 7.44 • 10“^AU. The excluded parameter region obtained 


from this experiment has already been derived in a previous work; see 


34( 1 for a full discussion. 


The second case we consider is that of a massless (or light) scalar field as discussed in section IVTl 
for which the PPN parameters approach the values (I82p . and are thus independent of the interaction 
distance. We may therefore also use bounds on the PPN parameters from experiments for which an 
interaction distance cannot be easily defined, such as the latest ephemeris releases INPOP13 |i 
331]. The bounds obtained from these datasets are given by 7 — 1 = (—0.3 ± 2.5) • 10“® and 
j3 — 1 = (0.2 ± 2.5) • 10“®. However, it turns out that the Cassini bound on 7 , and thus on a;, 
is still more stringent, and yields oj > 4.0 • 10^ at 2(7 confidence level. From the bound on /3 we 
then obtain the most stringent bound —2.5 • 10 < T — 4t(jcj < 2.7 • 10^*^ at 2cr confidence level for 

a; = 4.0 • 10^, and less stringent bounds for larger values of w. 

This concludes our discussion of the post-Newtonian limit for the general class of Horndeski’s 
gravity theories. The result we obtained can now easily be applied to particular theories within 
this class. We will show three examples in the following section. 


VIII. EXAMPLES 

After discussing the post-Newtonian limit of the most general form of Horndeski’s theory com¬ 
patible with our assumptions in the previous sections, we now come to particular example theories 
which fall into this class. In the following we will derive only the relevant PPN parameters 7 and (3 
for these theories, as we argued in the preceding section. In particular, we will discuss a common 
class of scalar-tensor theories with arbitrary potential in section IVHI Ai generalized Higgs inflation 
in section IVHI HI and the galileon model in section IVHI Cl 

A. Scalar-tensor gravity with a general potential 


As we already mentioned in the introductory section [U the work presented in this article gener¬ 
alizes a result obtained for a more restricted class of scalar-tensor theories 
of this earlier work is an action of the form 


34 ( 1 . The starting point 




dp(l)d^(j) - 2 k^V{(I)) + Sm[gi_LU,Xm], 


(84) 
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where the potential V{(j)) and the kinetic coupling function uj{4>) are free functions of the scalar 
field. Note that this action is written in the so-called Jordan frame, which is most convenient 
for calculating the PPN parameters and which can directly be compared to the action m we 
use in this article. Note, however, that one can also write this action in the so-called Einstein 
frame by application of a conformal transformation, and that the PPN parameters can also be 
calculated in this frame 35|]. This invariance of the theory under conformal transformations also 


allows expressing its PPN parameters in terms of invariants under these transformations 


361]. 


Since we use the action in the Jordan frame, we can directly compare it to the Horndeski gravity 
action ([1]) and read off the functions 


K{cP,X) = '^X-V{<P), G4(</.,X) = 


G3icP,X)=G5icP,X)=0. 


(85) 


K?(j) ’ ’ 2k‘^ 

We then expand these functions in a Taylor series around the cosmological background value cj) = ^. 
The relevant, non-vanishing coefficients in this Taylor series are given by 


^{ 0 , 0 ) = -^0 , ^(1,0) = -^1 ) ^(2,0) = -^2 , -fi'(3,0) = -^3 , 

Wo „ — Wo „ „ 1 


^(0,1) - 


^(1,1) = 


G, 


G, 


( 86 ) 

which takes 


^2^2 ’ 2^2 ’ 2k2 

Here we have expanded the functions w(i^) and V{4>) in analogy to the expansion ([S 
the form 

w((/)) =wo + wiV’ + C)(V^2)^ y(^) = yp + + ^ 2^2 ^^ ^ (37) 

where (p = ^ + ip. The constraint (I27p . which ensures the validity of the post-Newtonian approx¬ 
imation and its consistency with the cosmological background, takes the form Vq = W = 0, while 
the constraint (1501) is satisfied identically. For the constants determining the PPN parameters 
defined in equations (l33]) . (|i3]l and (fM)) we then find the values 


= 2k 


V2<^ 


o; = Wo , r = wi4> —Wo, a = 0 , /r = — 2 k^$^V3 . 


( 88 ) 


2too -|- 3 

Finally, we make use of these values in order to obtain the PPN parameters. These are given by 

2UJ0 + 3- e-™’/-'- 


j(r) = 


2uJo+ 3 + 


(89) 


and 


/?(r) = 1 + 


‘hwi 


(2^0 + 3-0 e-™v-^)2 I 2a;o + 3 


2m^r 


+ (2a;o + 3)m.^r 


e ln(m,/,r) — {m^r + Ei{—2m^r) — -e 


—2ra-^jJr 


+ x 1- 


$+3 


+ 


+2 2a;o + 3 


[e™'’^^Ei(—3m,;,r’) — e ™’^^Ei(—m,;,r)] 


(90) 
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which reproduces our previously found result. Finally, the massless case is given by V 2 = V 3 = 0 
and yields the PPN parameters 

Wo +1 a 1 , diwi 


7 = 


/3 = 1 + 


wq + 2 (2wo + 3)(2wo + 4)^ _ 

which is a well-known result for the PPN parameters of a massless scalar-tensor theory |49l | . 


(91) 


B. Generalized Higgs inflation 


We now discuss a class of models whose basic idea is to identify the Higgs field with the inflaton, 
which is the scalar field responsible for the inflation in the early universe, and which can thus be 
summarized under the name Higgs inflation models. It has been shown that a number of these 
models can nicely be written as a particular subclass of Horndeski’s gravity theory, which has been 


called generalized Higgs inflation 


5ril |. The functions in the action ([T|) of this model take the form 


X) = /C(0)X - V{ct>) , G3(./>, W) = h3(<())W , 

G4(<(-, X) = g{ct>) + h^X , G 5 (<(>, X) = h^{ct>)X 


(92) 


with six free functions JC, V, h^, h^, of the scalar field </>. After performing a Taylor expansion 
of these functions in analogy to the expansion ([8]) we find the relevant, non-vanishing Taylor 
coefficients 

^(0,0) = -bb , ^(1,0) = -^1 ) -^(2,0) = -^2 , ^(3,0) = -^3 , -^(0,1) = '^0 , ^(1,1) = ) 
G^3(0,1) = ^3,0 ) ^4(0,0) = 90 ) G^4(1,0) = 9l ) ^4(2,0) = 92 , ^4(0,1) = ^4,0 ; ^^4(1,1) = ^4,1 • 

(93) 

Similarly to the scalar-tensor theory discussed in the previous section we find that the con¬ 
straint (EH), which ensures the validity of our perturbative expansion around a stationary cos¬ 
mological background, takes the form Vq = Pi = 0. Further, the constraint (isni) translates to 
^3,0 = 3/i4^i and = 0, which we will impose in the remainder of this section. We then find the 
mass of the scalar field and the constants 

_ 9o92 _ 9 q^i 


= 


‘^90^2 


UJ = 


a = 


9o92 
2 ’ 
91 


T = 




3 ’ 


= 


9SVs 


9i 


(94) 


y 90^0 + Sgf ’ 2gl 

from which the PPN parameters 7 and B can be obtained. 

□ 

We discuss a few special cases listed in [^, which can be viewed as corrections to general 
relativity with a minimally coupled scalar field given by the gravitational action 

/M2 


Sg = I + X- ViB) + 


(95) 
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where Mpi is the Planck mass. Without any such correction A/i the only non-vanishing Taylor 
coefficients are go = Mpj/ 2 , /Cq = 1 and the terms originating from the potential In 

particular, we consider the following models: 


Running kinetic inflation [^, AC = with parameters k and n. In this model we 

obtain the modified Taylor coefficients /Cq = 1 + and KLi = In this case we 

have gi = 0, so that we obtain the limit w —>■ oo, from which follows 7 = /3 = 1. 


Higgs G-inflation AC = —(pXdcp/M'^ with parameter M. In this model the only 

modified Taylor coefficient is However, this coefficient is restricted by the 

condition (|5np . so that this model will yield terms in the gravitational field equations which 
are not covered by the PPN formalism we used in this article. 


Non-minimal Higgs inflation: 


5445611 : AC = —^(lP‘R/2 with parameter In this model we 


find the modified Taylor coefficients 


90 = 


Mil - 


91 = , 92 = - 


e 

2 ’ 


(96) 


from which follow the constants 


Tfl^ — Meff .4 


2 P 2 


M2j + 6C2^>2 ’ 


U) = 


^eS 

4 ^ 2(^2 


a = — 


^eS 

4^4.2 


r = 0 , /i = 


4^3 (^3 


(97) 


where we introduced the effective Planck mass M|j = M|i — Typically one is interested 
in the case Mgfj ~ Mpj. Here we make use of the fact that Mpi 3> ~ 246GeV, which is 

the vacuum expectation value of the Higgs field, and expand the PPN parameters 7 and /3 
in orders of ‘h/Mgff. Up to the first non-trivial order we find 


7 


= 1 - 


(J)2 


+ 0 


$3 

w« 


(98) 


and 


/3 = 1 + ^- 26 “”^’^^ In(m^r) 


+ 2(m^r + e-^'')Ei(-2m^r)] + O . (99) 


New Higgs inflation 57H6fll| : 


AT = ^ [XR + (00)2 - 


( 100 ) 
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with parameter jl. In this model the only modified Taylor coefficient is /i 4 ^o = 1/2/i^- 
However, as it was also the case for Higgs G-inflation, this coefficient is restricted by the 
condition (1501) . so that also this model will yield terms in the gravitational field equations 
which are not covered by the PPN formalism we used in this article. 


• Running Einstein inflation 



A 6 


1 ^- ^3 , w .^2 1 


- -{a<pr + -(□0)(v^v,0)^ - 


( 101 ) 


with parameter A. This term does not influence any relevant Taylor coefficients, so that 
again we find 7 = /3 = 1 . 


We thus see that the non-minimal Higgs inflation model is the only model from which we obtain 
PPN parameters which potentially deviate from observations. Note, however, that for a scenario 
in which ~ 125GeV is the Higgs mass, any observable gravitational interaction takes place in 
the limit m^r S> 1 , such that one has the limiting values 7 = /3 = 1 . 


C. Galileons 


The last example we consider here is a scalar field whose action, in flat spacetime, is invariant 
under Galilean transformations, and hence is called galileon 
theory in curved spacetime 6 ^. The action is given by 


6 l|. Here we consider the covariant 


S= d^x^/^\R + C4Vx4>^^ 


1 . 


2incj)y - 2(V^V,,/.)" - 


+ -c5Va</>vv [(□</>)" - + 2(v^v,</.)’ 


( 102 ) 




After integration by parts one can see that this action has the form of the Horndeski action p!|). 
where the free functions are given by 


K{^,X) = Ci^-2C2X, G3{cP,X) = 2C3X, ^ 4 ( 0 , X) = 1 - 2 C 4 X 2 , G5i^,X) 

The relevant, non-vanishing Taylor coefficients are thus given by 


ISCsX^ . 

(103) 


-^( 0 , 0 ) — ) -^( 1 , 0 ) — C”! , K(^q i'j—— 2 C 2 , ^ 3 ( 07 ) — 2 C 3 , ^ 4 ( 0 , 0 ) — 1 • 


(104) 
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We see that the terms involving (74 and do not enter the relevant Taylor coefficients, and so 
have no influence on the post-Newtonian limit. The conditions (1271) and (1501) . which we imposed in 
order for the post-Newtonian limit to be valid, require that (7i = (73 = 0. From 77(2,0) = -^( 3 , 0 ) = 0 
further follows that the scalar field is massless and we can apply the formalism detailed in section IVTl 
It turns out that 7 = /3 = 1 for this class of theories. 

This concludes our discussion of particular examples for Horndeski gravity theories. We have 
seen that among this class there are several theories whose post-Newtonian limit is consistent with 
the solar system observations displayed in section IVIIl as they are in particular compatible with 
the values 7 = /3 = 1 obtained for various example theories. 


IX. CONCLUSION 

In this article we discussed the post-Newtonian limit of Horndeski’s theory of gravity. We 
showed that the post-Newtonian limit is fully determined by fifteen constant parameters, which 
arise as coefficients in the Taylor expansion of the free functions 77, ^ 3 ,^ 4 , G 5 in the Horndeski 
Lagrangian around the cosmological background value of the scalar field. It turned out that 
for the post-Newtonian limit to be consistent, we must impose several constraints on these Taylor 
coefficients. With these constraints in place, we calculated the post-Newtonian limit in two different 
scenarios. 

In the first scenario we considered the most general theory consistent with the aforementioned 
constraints. We defined the PPN parameters 7 and /3 and calculated their values for a static 
point mass. It turned out that the PPN parameters are not constant, but depend on the distance 
between the gravitating source and the test mass. 

In the second scenario we imposed additional constraints on the Taylor coefficients under which 
the scalar field becomes massless. We showed that when these constraints are satisfied, the post- 
Newtonian limit of the theory assumes the standard PPN form, which is characterized by ten 
constant PPN parameters. We calculated these parameters and showed that Horndeski’s theory 
is a fully conservative theory, which is free of preferred-frame and preferred-location effects, which 
means that only the parameters 7 and /3 potentially deviate from their observed values. 

We finally applied our analysis to a number of example theories, including a previously discussed 
scalar-tensor theory with arbitrary scalar potential, various Higgs inflation models and galileons. 

The work presented here allows for further extensions and generalizations. A straightforward 
generalization is to drop the assumption that the cosmological background value of the scalar 
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field is constant and to allow for a time dependence <1> 7 ^ 0. Another possibility is to investigate the 
parameterized post-Newtonian limit of more general scalar-tensor theories beyond the Horndeski 
Lagrangian, which introduce higher order derivatives into the gravitational field equations. Despite 
originally being regarded as ill-defined due to Ostrogradski instabilities and ghosts, it has turned 


out that these problems may be overcome and healthy theories exist 


63h 69I|. Yet another possible 


direction of future research is to consider theories with several scalar degrees of freedom, such as the 


recently developed generalization of Horndeski’s theory to two scalar fields 


Q- 


Finally, one may 


also consider modifications of the formalism itself, in order to include effects caused by screening 
mechanisms such as the Vainshtein mechanism [7l|. 


Appendix A: PPN expansion coefficients 

This appendix lists the coefficients appearing in several lengthy equations and intermediate 
results of the calculation of the PPN parameters in sections and IVIl The following coefficients 
appear in the scalar field equation at the fourth velocity order: 


do = 


do = - 


2iF('o n'lG*. 


(-^(1,1) “ 4^3(2,0))G^4(0,0) + (-^(0,1) “ 2G3(i^0) + 12^4(2^0))G'4(1,0) 


2(.f^(0,l) “ 2G3 (i^0))G'4(0,0) + 6^4(1,0 ) 


(2,0)'-^4(1,0) 


— 3K(o n'l G, 


(3,0) <-^4(0,0) 


(-f^( 0 ,l) “ 2G3 (i^o))G'4(0,0) + ’ 


di = — 


(-^(1,1) “ 4^3(2,o))G4(0,0) + 6^4(1,0) 1^4(2,0) 


(-^( 0 , 1 ) “ 2(73(1^0) )G4(o,0) + 3^4(1, 0 ) 


d^ — —df\ — 


( 1 ^ 3 ( 0 ,1) “ 3^4(1^!) )G4(o,o) — (<^4(0,1) “ <-^5(1,0) )G4(i,0) 


ds — 


*^4(1,0) “ <^4(0,0) <^4(2,0) 


(-^( 0 , 1 ) — 2G3(i^0))G^4(0,0) + 3^4(1 0 ) 

2(04(0,1) “ 05 ( 1 ^ 0 ))O4{0,0) 


(-^(0,1) “ 2G3(i^o))04(o,o) + 30^1.^^ ^^ 


dm — 


di2 — —3dn — 


3G, 


4(1,0) 


2 (Ar(o,i) — 203 ( 1 ^ 0 ) )04(o,o) + 


(-^( 0 , 1 ) “ 2G3(i^0))O4(0,0) + 

di = dy = dg = 1 . 


(Al) 


From this equation one derives equation ()5ip . which contains the following coefficients: 


62 = 


327r2 


“Ido + (dy -|- d8)ci — (dy -|- 3d8 — 2dg)c^ -|- 2(d3 -|- d4)ni^ -|- 2dsni^ 


2 04 ( 1 , 0 ) 

04 ( 0 , 0 ) 


63 = esm^ = 3^2^^ ^7 [(c 3 - ci)c^ -h (c2 - C4)m^] , 


ei = dg (c4m^ - C3C^) , 64 = eem^ = [dy(ci - C3) -b 2d3m^] 


M^c 


(A2) 
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In terms of these coefficients we express the coefficients of the solution (|52p in the form 


“ 2m^ "^2 4 ’ 

. _ 65 j _ 66 . _ 

ji 2 ’ 2 ’ 


_ 62 _ ^ Oeem^ 

“ 2m^ 2 4 


6l 


2m. 


ip 


fb — 


A similar expansion is used for the metric component . The coefficients in equation 
listed below: 


6l 


2m. 


"Ip 


63 + 65777 ,^ 


(A3) 


are 


<^4(1,0),^ ,, „ G'4(i^o)6^2 - A'( 2 , 0 ) 

Ql = 7^ - (,“1 - -Lj > q2 = - 


Qi = 


^4(0,0) 

^4(1,0)6^4 + 2^4(2, 0 ) 


a 


4(0,0) 


95 — 


^4(0,0) 
^4(1,0) 


, 93 = 

^4(1,0) 


G^4(1,0)6^3 + 2^4(2, 0 ) 


a 


4(0,0) 

G, 


1 - j *^4(1,0) / , IN 

^ «5 , 95 — TP -“6, g? — 77-(«7 - 1) 

<-^4(0,0) 6x4(0,0) 6x4(0,0) 


G , 


98 = 


4(1,0) 


G . 


4(0,0) 


ds , 99 — 


G . 


4(1,0) 


G . 


4(0,0) 


(dg — 1), 9io — 


G . 


4(1,0) 


G , 


4(0,0) 


6^10 5 9 ll — 0 , 


G . 


912 = 


4(1,0) 




913 = 2 


G 


4(1,0) 


”^4(0,0) ^4(0,0) 

2^4(1 0 ) dll ~ 1 1 

917 = - 777 -) 918 = 


gi4 = -l, 915 = 1, 916 = 


G, 


4(h0)_2 


m ,. 


ri ■■''Ip 1 
6^4(0,0) 


, 919 — 


1 2G4(i o)di2 ~ 3 

, 920 = 


2G4(o,o) G4(o,o) 2G4(o,o) 2G4(o,o) 

From this equation one obtains equation (|55l) . which contains the following coefficients: 


(A4) 


•52 = 916/2 + 


327r2 


2 q 2 + 2(^3 + 9 i)m‘i + ^7(03 - ci) + q's ci - 3 c 3 + 2 


G, 


4(1,0) 


G, 


Cip 


+ 299 C 3 + ( 2 g'ii + 2 g'i 2 — 913)61 + 2 


91461 + 915C1C3 


m 


'ip 


Si = 


M2 


914 62 


4(0,0) 


CiC^p 

m^p 


53 — S^Ulip — 


•Si = 916/1 + iq^ 2^2 + 9i5)(c47n^ “ 636^) + 9i2(c2m^ - Cic^)m^] , (A5) 

3^^2^ [97((c 3 - Ci)c^ + (C2 - C4)m^)m^ - 2(giim^ + 2gi4Ci)(cic^ - C2m^)] , 


54 — ^SQiri'ip — 


M 24 


IGvr^ 


m 


'ip 


[293m^ + qrici - 63)771^ + 2 giicim^ + 2^1464] , 


•S8 — 910/3 , Sg — 910/4 , Sio — 910/5 , Sll — 916/6 
The solution (|56p for is given in terms of the following coefficients: 


S 8 + Sio '55711^ Ss — SlO 

uj — S3 T s^TTipj 2 , Tig — Ui^mpj — Si S3Tnp) T — -t- 


m 


Ip 


mpj 


S9 “t“ <811 2 <89 — 811 , , 

uio = -84 + ‘2'SQmp, -—2— , till = 2 it 5 m^ = 82 - 2s4mp, + 2sem^ H-^3^-, (A6) 
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Ip 


mp, 


87 


85 


Y’ ^^- 2 ' 


S 6 S 8 59 510 

tt3 — — 5 ^6 — —2 ? ^9 — —2 ’ ^12 — —2 ’ ~ 


5ll 


mt;. 


mz 


m:, 
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In the case of a massless scalar field one obtains equation (17211 with the following coefficients: 

2gi9C2 - 2q^cl - qgc^ (^C2 - 3c4 + - 2qgc^c^ - {2qi2 - gi3)c2C^ - 2^150204 


tcs = 


Stt 

2q'34 + 97(c 2 - C4)c^ + 2giiC2C^ + 2gi4C^ 
32 ^ 1 ^ ^ 

Wi = -r^ , Wi = qn , rcs = g'ls , wq = ^20 
dvr 


W2 = 


(A7) 
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